Experimentally realizable characterizations of continuous variable Gaussian states 
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Measures of entanglement, fidelity and purity are basic yardsticks in quantum information pro- 
cessing. We propose how to implement these measures using linear devices and homodyne detectors 
for continuous variable Gaussian states. In particular, the test of entanglement becomes simple with 
some prior knowledge which is relevant to current experiments. 
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In the current development of quantum information 
processing, we have witnessed the importance of fidelity, 
entanglement and mixedness. In this paper, we propose 
feasible experimental schemes to measure these critical 
quantities for Gaussian continuous- variable systems. The 
proposed schemes require standard laboratory devices 
such as beam splitters and phase shifters and highly effi- 
cient homodyne detectors. 

Entanglement has been mainly confined to theoretical 
discussions. Only very recently, Horodecki and Ekert 
investigated how to measure entanglement. Experimen- 
tal studies on quantum correlation between two or more 
particles have been concentrated on tests of the Einstein- 
Podolsky- Rosen (EPR) paradox and Bell's inequality be- 
cause an experimental measure of entanglement was not 
so clear. The present paper proposes an experimental 
scheme to measure a degree of entanglement with some 
prior knowledge of the given state. There are some indi- 
rect ways to test entanglement, for example, by proving 
the fidelity higher than the classical limit in quantum 
teleportation p|. However, this provides only a suffi- 
cient condition for entanglement. On top of that, how 
to measure the fidelity has not been thoroughly investi- 
gated. Here, we also propose an experiment to measure 
how close two states are. 

Experiments on continuous-variable quantum infor- 
mation processing have been concentrated on Gaussian 
states ||, ||, |j. This is because, due to extremely low 
efficiency and high dispersion in high-order nonlinear in- 
teraction, experiments have been based on linear trans- 
formation of fields initially in thermal equilibrium. In 
thermal equilibrium, a field is in a Gaussian thermal 
state H and the linear transformation keeps the Gaus- 
sian nature ||. The linear transformation of an input 
field to an output fields is due to a Hamiltonian com- 
posed of linear and/or quadratic bosonic operators a and 
a*. Any linear transformation for two-mode fields can 
thus be represented by a product of single-mode squeez- 
ing, rotation and displacement operations and two-mode 
squeezing and beam-splitting operations 0]. 

To measure how close a quantum state p\ is to a ref- 
erence pure state |f/>) 2 , the fidelity is defined as T = 
2{\pi\ip)2- The fidelity is important, for example, to find 
how successfully a state is reproduced after a set of local 



quantum operations and classical communications such 
as a teleportation process. However this theoretical con- 
cept has not been thoroughly investigated by experiment. 
We propose a feasible experimental scheme to realize the 
measurement of the fidelity. 

It is convenient to work with the Weyl characteristic 
function defined as Ci(x) = Trpii)j(x). Here, the dis- 
placement operator is defined as -Di(x) = exp(ixiX T ), 
where the operator vector Xj = (qi,Pi) with quadrature 
operators <j, = (hi + dj) and pi = i(a\ — hi) and the co- 
ordinate vector x = (p, —q). The fidelity T is equivalent 
to the overlap between the characteristic functions: 



T 



1 



dxCi(x)C 2 (-x). 



(1) 



Throughout the paper, matrices are represented in bold 
face and operators with hats. In order to measure the 
fidelity of two fields, they are mixed at a beam splitter 
whose action is described by Byi — exp[#(djd2 — didj)], 
where the reflectivity r and transitivity t of the beam 
splitter are determined by 8: r = sind and t — cos 8. 
The characteristic function for an output field is 

C out (x) = Tt 2 £ 2 (x)'&i£i2ftA! - Ci(rx)C 2 (tx). (2) 

By Fourier transforming the characteristic function, the 
Wigner function, W(a), is obtained M. Thus we can 
easily see that the Wigner function of the output field at 
the origin (a = 0) of phase space is directly related to 
the fidelity: 



T= — \ dxCovtW =2ttVF(0) 
Ztt 



(3) 



when r = — t = l/v2- We have found that, after mixing 
the two fields at a 50:50 beam splitter, we measure the 
Wigner function at the origin of the phase space for one 
of the output fields, to find how close the two input fields 
are. The Wigner function of a given field can be mea- 
sured using optical tomography [[7|, which requires some 
numerical processes on experimental data. However, if 
both the input fields are Gaussian, as the output field 
is also Gaussian, W(0) is easily measured using a highly 
efficient homodyne detector. 
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The characteristic function of any single-mode Gaus- 
sian field is written as 

C(x) = exp(-ixV s x T + idx T ) (4) 

where V s is the 2x2 variance matrix defined as (V s )ij = 
h(({xi, Xj})) — (xi)(xj) and d is the displacement vector 
di = (xi). The quadrature variables are measured by a 
balanced homodyne detector B, which is a well-known 
device to detect phase-dependent properties of an opti- 
cal field. The operational representation of the balanced 
homodyne detector is Ohd = qcos\ + ips^Xi where % 
depends on the local-oscillator phase. The value W(0) of 
the Wigner function is rotationally invariant so the nu- 
merical process becomes simpler as we rotate the output 
field to diagonalize the variance matrix V s . This can be 
done by placing an optical phase shifter before the homo- 
dyne detector. When the variance matrix is diagonalized 
the uncertainty AqAp is minimized and the fidelity (|^) 
is 



f 



AqAp 



exp 



2\Aq 2 Ap 2 



(5) 



where (p), (q), (p 2 ) and (q 2 ) are measured by a homo- 
dyne detector. We have shown how the fidelity for two 
Gaussian fields is measured using a homodyne detector 
and a beam splitter. In relation to the two input fields, 
(p) and (q) are proportional to the displacement differ- 
ence between them and AqAp is the uncertainty over 
their averaged variance matrix. When the displacement 
difference for the two fields is zero, the fidelity is deter- 
mined by their average uncertainty: T = 1/ApAq. 

We have proposed an experimental scheme to measure 
the fidelity of two single- mode fields. We are now inter- 
ested in how to measure the a degree of entanglement 
for a two-mode Gaussian field. In the discussion, we 
will be able to show how the entanglement measure is 
related to the mixedness of the system. The entangle- 
ment of a two-mode state does not change by displace- 
ment operation j^] so we write the general form of the 
characteristic function, C(x), for a two-mode Gaussian 
state without the linear displacement term in Eq. (|4|): 
C(x) = cxp(— xVx T /2), where x is now ?2>P2) 
and V is a 4 x 4 variance matrix for modes 1 and 2. The 
variance matrix can in fact be written using 2 block ma- 
trices Li and L2 for local quadrature variables and C and 
its transpose C T representing inter-mode correlation, 



V = 



Li C 
C T L 2 



(6) 



What are the possible types of entangled continuous- 
variable states one can produce? We start with two in- 
dependent fields of modes 1 and 2, which are in ther- 
mal equilibrium at temperatures T\ and T 2 , respectively. 
The variance matrix has only local elements, Li = nil 
and L2 = U2lj where 1 is the 2x2 unit matrix and 
hi = 1 + 2/[exp(huj/kBTi) — 1] with the Boltzmann 



constant ks- The basic components of a linear trans- 
formation are single-mode squeezing, Sj(s), and rota- 
tion, Rj(0), and two-mode squeezing, Si2(s), and beam- 
splitting, B12, whose actions are described by their trans- 
formations of the variance matrix: 
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(7) 



where er z is a Pauli spin matrix. Any combination, U, 
of these matrices, transforms the variance matrix V into 
U T VU. Here, displacement has not been considered be- 
cause it does not affect entanglement. 

The degree, s, of squeezing normally determines the 
optimum entanglement for a given setup [^). As squeez- 
ing is due to x^ 2 ) -nonlinear interaction, the efficiency 
is relatively low and a combination of squeezers is not 
practical |I(J. Entangled states have been produced us- 
ing two- mode squeezing S12 |3| or the combination of 
Bi2Si(r)S2(— r) 0. For these linear operations, the vari- 
ance matrix is represented by diagonal Li, L2 and C. In 
particular, if the two modes are initially in thermal equi- 
librium at the same temperature before a transformation, 
the local matrices Li and L2 are identical, for which case 
we show that the degree of entanglement becomes simple 
and can be measured using joint homodyne detectors. If 
such a Gaussian field is decohered in two-mode thermal 
bath of which the two modes are in the same tempera- 
ture, the decohered state is still represented by the same 
form of variance matrix . 

The inseparability of a continuous variable Gaussian 
state is determined by Peres- Horodecki condition J12| . A 
density operator p is entangled if and only if the partially 
transposed density operator fF 2 has any negative eigen- 
values. In order to quantify entanglement, we define the 
degree of entanglement, £(p), as the absolute sum of the 
negative eigenvalues of the partially transposed density 
operator: £(p) = Tr|/5 T2 | — 1. This is an entanglement 
monotone |13[ . In the following, we show how £(p) is 
related to homodyne measurements. 

Lemma 1. - If the block matrices Li = L2 and C 
are diagonal, i.e., the variance matrix of a Gaussian 
continuous-variable state has the following form: 



V 



fm a 

n 2 c 2 

ci Ui 

V c 2 n 2 j 



(8) 



where ri\ or n 2 may be smaller than the vacuum limit 1. 
The degree of entanglement is given by 



£{p) = max{0, (S^y 1 - 1} 



(9) 



where Si — n, — |cf| for i = 1, 2. 

We briefly sketch the proof. The main task is to cal- 
culate the trace of |/5 T2 | which is the positive operator 
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satisfying \p 



T 2 |2 



(p 2 ) . In order to find the bound 



operator \p 2 \ and its trace, all operators are represented 
by their characteristic functions respectively, based on 
the one-to-one correspondence principle between a bound 
operator and its characteristic function Q . Then the op- 
erator equation, \p T2 \ 2 — (p T2 ) 2 , is converted into an 
equation for their corresponding characteristic functions. 
The equality of the two Gaussian characteristic functions 
implies that their variance matrices and normalization 
values are the same. Let V p and N p be the variance 
matrix and the normalization value for \p T2 \. The unit 
trace of p T2 leads its normalization value to be unity: 
N = Trp T ' 2 = 1. Based on the observation that the 
transposition is momentum reversal, the variance matrix 
Vo of fF 2 is obtained from Vo with p 2 — > —p2 ■ One thus 
has the two equations for V p and N p as 



nv p 1 n T = v - raVo L n 



V detV 



(10) 



(11) 



In addition, because \p T2 \ is positive, its variance matrix 
V p satisfies the uncertainty principle ]l2] | , 



V B + ifl > 0. 



(12) 



where l~i = icr,, 



1 with the Pauli spin matrix oy From 

in 



the solution V p satisfying both Eqs. (|10| ) and (|12|), N p i 
Eq. ([fl]) is obtained. Finally, the degree of entanglement 
£ = Nj> — 1 as Tr\p T2 \ — N p and £ is given in particular by 
Eq. (m for Vo. If and only if the partial transposed den- 
sity operator is positive satisfying the uncertainty princi- 
ple ([l2]), V p = Vo and the state is separable with £ = 0. 

For the two-mode vacuum, n± = n 2 = 1 and c± = 
C2 = so that (5i)o = 1 where (-)o denotes the value 
for the vacuum. Recalling the definition of the vari- 
ance matrix, <5i = \((q\) + (gf) - 2|(g 1 g 2 )|) and 5 2 = 
\((p 2 ) + (pi) — 2|(pip 2 )|), which can be measured by 
joint homodyne detectors for modes 1 and 2. Accord- 
ing to Eq. (||) the state is separable with £ = when 
S1S2 > {6162)0- Otherwise, the degree of entanglement 
is £(p) = (6162)0/(6162) — 1- We have found that the 
entanglement of a Gaussian field in the form (|J) can be 
tested by comparing the joint quadrature variance of the 
given field with that of the vacuum. 

Let us consider the purity of a state. A degree of purity 
can be defined as P = Trp 2 = (1/2tt) / dxC(x)C(-x), 
where C(x) is the characteristic function of the given 
state. When P = 1 the state is pure. The purity of a 
Gaussian state with its variance matrix Vo is easily cal- 
culated using the determinant of Vo: P = l/\/dct Vo. 
We have already seen that the matrix elements of Vo 
can be measured by homodyne detectors so the purity of 
the Gaussian state can be measured. To find the rela- 
tion between the purity and entanglement, we define the 
mixedness M of a system: M = 1 — P, which is zero when 
the system is pure and becomes unity when it is totally 
mixed. According to Lemma 1, a Gaussian state of the 
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FIG. 1: Configuration of the setup to test entanglement. The 
dotted devices measure off-diagonal terms of the local vari- 
ance matrices Li and L2. The boxed devices denoted by HD 
are homodyne detectors. The numbers refer to the modes. 



variance matrix Vo is separable when 6162 — 1 > 0. Mul- 
tiplying a positive quantity (m + 1 ci | ) (n 2 + 1 c 2 1 ) — 1 to the 
both sides of this inequality, we find the necessary and 
sufficient condition for separability: M12 — Mi — M 2 > 
2|cic 2 |, where Mi and M 2 are the mixedness of the field 
in mode 1 and 2 and Mi 2 is that of the whole field. For a 
pure two-mode state with Mi 2 = 0, the state is separable 
iff Mi = and M 2 = with \ac 2 \ = 0. 

So far we have assumed that the block matrices Li 
and L 2 and C are diagonal because this case is relevant 
to the current experimental conditions and the degree of 
entanglement becomes extremely simple. However, how 
do we make sure that there is vanishing off-diagonal ma- 
trix elements? Measuring the off-diagonal terms of local 
matrices is troublesome as it involves the joint measure- 
ment of two quadrature variables for a single mode. To 
measure the off-diagonal elements of Li Q, we put a 
50:50 beam splitter which splits the field in mode 1 as 
schematically shown in Fig. fil. Using the beam splitter 
transformation in (^) with r = t = 1/ \/2, we find that 
the field for three modes 1', 2 and 3 is still Gaussian and 
its variance matrix is written as 



L! + l V2C -Li + lP 
V2C T L 2 -V2C T 
-Li + 1 -V2C Li + 1 



(13) 



where the unit matrix 1 is due to the vacuum injected 
into the unused port of the beam splitter. Now the off- 
diagonal elements of Li can be measured by inter-mode 
correlation between modes 1' and 3: The mean value of 
the joint measurement (qvP3) = — V12/2. Similarly other 
off-diagonal terms of the local variance matrices can be 
obtained. 

In fact, we have shown how to find all the matrix ele- 
ments of the variance matrix for a Gaussian field so that 
it is possible to test entanglement not only for the fields 
in the form (|J) but also for any Gaussian field if the 
detection efficiency is unity. In this case we need to gen- 
eralize the expression (^) for the degree of entanglement 
which is rather straightforward using the same argument 
to derive (^). An equivalent but alternative approach 
was suggested by Vidal and Werner [jl3| . 

A homodyne detector is composed of two photodetec- 
tors. Inefficient photodetectors may miss photons to de- 
tect and reduces the quantum correlation between two 



4 



modes. The detection efficiency thus determines the fea- 
sibility of the proposed schemes, in particular, to test 
entanglement of the field. If the efficiencies of the pho- 
todetectors are same, homodyne measurement by imper- 
fect detectors is equivalent to homodyne measurement 
by perfect detectors following a beam splitter, one input 
port of which is fed b y t he field to be measured and the 
other by the vacuum |lq| . Here, we assume non-unit de- 
tection efficiency due only to missing photons to detect. 
The efficiency r] of the homodyne detector determines 
the transmission coefficient ^fr\ of the beam splitter. In 
fact the fictitious beam splitter affects the testing field as 
though it is decohered in the vacuum reservoir. The de- 
tection efficiency, assumed the same for the both homo- 
dyne detectors, effectively changes the variance matrix 
from V to V = 77V + (1 — This is what is measured 
by imperfect homodyne detectors. If V is in the form 
Vo, the variance matrix Vq takes the same form as Vo 
but with modified matrix elements n! i = ryrii + (1 — rf) and 

\ c i\ = v\ c i\ f° r eacn i- 

Consider the effect of the detection efficiency on the 
inseparability of the testing fields. Substituting n! i and 

into the separability condition, > 1, for inefficient 
detection, we find a state to be entangled when 

(»' - V1W2 - 141) = ?? 2 (<W2 -5,-52 + 1) 

+r)(S 1 +S 2 -2) + l<l. (14) 

Rearranging this equation, we can easily find that when 
the original testing field is characterized by 5, + 62 < 
2, it is always found to be entangled regardless of the 
detection efficiency unless the efficiency is zero. Fig. || 
presents the sets of Gaussian states on the space of 5\ 
and 62 where separable states with 5\52 > 1 are denoted 
by S and entangled states with 5i5 2 < 1 by E + E' . 
All entangled states in the region E with the condition, 
§1+62 < 2, will violate the inequality, 5\52 > 1, unless the 
detection efficiency is zero while some entangled states in 
E' fail the test of entanglement. 

Reid and Drummond derived the inequality for the 
quantum correlation between two mode fields along the 
line with the EPR argument ||, [l(|. They introduced 
the uncertainty V\ (V 2 ) between the observable q\ (pi) 



in one mode and q[ (p^) inferred from the observation 
of the other mode. Quantum correlation may lead the 
product of the uncertainties to be less than the vacuum 
limit, resulting in the inequality of V1V2 < 1. In our 
notation this inequality can be written as 

x x ^ 711712 

(n x + | Cl |)(n 2 + |c 2 |) 

Note that the right hand side of the inequality is always 
less than unity. Thus, the inequality (|l5| ) is sufficient to 
satisfy our inseparable condition (||) JlT]]. However, the 
converse statement does not hold in general. 

We have proposed experimental schemes to measure 
the fidelity, the purity and the degree of entanglement of 
a given system with some prior knowledge. The scheme 
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FIG. 2: Gaussian states with the symmetric variance matrix 
Vo in Eq. (|^). Separable states are denoted by S and en- 
tangled states by E and E' with the boundary of 8182 = 1. 
Further the two regions E and E' are separated by the line 
<5i+<52 = 2. The entanglement imposed on a state of the region 
E' may fail the entanglement test due to inefficient detection. 
However, the entangled state of the region E is always found 
to be entangled regardless of the detection efficiency. 

consists of beam splitters and balanced homodyne de- 
tectors which are well-established experimental tools to 
study quantum optics. When the detection efficiency is 
unity, the scheme can be extended to a general two-mode 
Gaussian state to test its entanglement. 
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